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Detailed  kinetic  calculations  of  chemical  reactions  using  GRI-Mech  3.0  are  carried  out  to  characterize  the 
combustion  of  CH4/02/N2  mixture  in  a  well-stirred  reactor  (WSR).  Very  wide  ranges  of  the  inlet  temperature 
Tin  (300  K-2700  K),  nitrogen  dilution  concentration  (0%— 99%)  and  global  equivalence  ratio  (0.5-4.0)  are 
considered.  The  extinction  and  self-ignition  temperatures  (Tex  and  TSi)  of  the  CH4/02/N2  mixture  are  iden¬ 
tified  by  the  ignition-extinction  S-curve.  Based  on  lex*  Ts  :i  Twsr.  and  the  operative  conditions  (Tjn  and  the 
mixture  composition),  the  WSR  combustion  of  CH4/02/N2  can  be  quantitatively  classified  into  several 
particular  regimes.  Moreover,  the  product  composition  and  elementary  chemical  pathways  of  the  CH4 
oxidation  are  also  examined.  Results  demonstrate  that  the  WSR  working  temperature  determines  the 
chemical  characteristics  of  the  CH4  oxidation  rather  than  the  combustion  regime. 
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1.  Introduction 

The  moderate  and  intense  low-oxygen  dilution  (MILD)  com¬ 
bustion,  termed  formally  by  Cavaliere  and  de  Joannon  1  ],  is  one  of 
the  newly  developing  combustion  technologies.  This  combustion  is 
similar,  but  not  identical  (see  below),  to  those  that  debuted  a  few 
years  earlier,  i.e.,  the  flameless  oxidation  (FLOX)  [2  and  the  high 
temperature  air  combustion  (HiTAC)  [3,4  .  So  far,  significant  ben¬ 
efits  of  these  technologies  for  industrial  applications  have  been  well 
demonstrated,  e.g.,  producing  efficient  combustion  4,5)  and  in 
particular  extremely  low  NOx  emissions  [6,7].  These  technologies 
are  commonly  featured  by  highly  distributed  reaction  zone  [8], 
uniform  thermal  field  [9],  and  low  luminescence  [4,10  .  Accord¬ 
ingly,  in  recent  years  the  HiTAC  and  FLOX  are  often  treated  identical 
to  the  “MILD  combustion”.  Note,  however,  that  those  technologies 
are  not  exactly  identical,  as  indicated  by  Cavaliere  and  de  Joannon 
[1  and  further  argued  below. 

Wiinningand  Wiinning  2]  carried  out  a  series  of  experiments  on 
diffusion  combustion  of  burning  methane  (CH4)  using  a  flameless 
oxidation  (FLOX)  burner.  Based  on  the  data  obtained,  they  delivered 
a  mapping  of  different  combustion  regimes.  As  reproduced  in 
Fig.  1(a),  three  combustion  regimes  are  evident:  (i)  traditional 
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combustion  (TC)  occurring  possibly  over  the  entire  range  of  the 
furnace  temperature  but  only  for  the  dimensionless  recirculation 
rate  I<v  <  0.3  (somewhat  increased  at  higher  temperatures);  (ii) 
flameless  oxidation  (FLOX),  i.e.,  combustion  reactions  occurring 
invisibly  in  a  highly  stable  form  when  both  the  furnace  temperature 
is  well  beyond  the  self-ignition  temperature  of  the  mixture  (Ts 0  and 
the  entrained  flue  gas  recirculation  rate  Kv  >  3;  (iii)  unstable  com¬ 
bustion  (UC)  for  3  >  Kv  >  0.3,  which  intermittently  switches  be¬ 
tween  TC  and  FLOX.  Within  the  UC  regime,  the  flame  may  lift  off  and 
finally  extinct  for  temperatures  <  Ts + 

Nippon  Furnace  Kogyo  Kaisha  Company  3]  meanwhile,  taking  a 
different  route,  developed  a  new  combustion  process  by  preheating 
the  combustion  air  to  a  temperature  Ta[r  >  Ts j,  usually  around 
1200  K.  Katsuki  and  Hasegawa  [3]  termed  it  as  the  high  tempera¬ 
ture  air  combustion  (HiTAC),  where  the  “high  temperature”  means 
that  rair  >  TS[.  Namely,  HiTAC  is  somewhat  different  from  the  FLOX, 
see  Fig.  l(a,  b).  Tsuji  et  al.  [4]  considered  the  traditional  preheated 
air  combustion  (PAC)  as  the  one  with  Tajr  <  Ts j.  Fig.  1(b)  shows  the 
auto-ignition  and  forced-ignition  limits,  derived  from  the  oxygen 
concentration  and  air  temperature  (two  external  injection  param¬ 
eters),  for  propane  (C3H8)  jet  in  hot  air  diluted  with  nitrogen,  from 
which  two  combustion  regimes,  i.e.,  PAC  and  HiTAC,  are  defined. 
Comparison  of  Fig.  1(a)  and  (b)  reveals  that  the  traditional  and 
unstable  combustion  (TC  and  UC)  at  Tair  >  TS[  and  FLOX  of  Wiinning 
and  Wiinning  2]  altogether  are  equivalent  to  the  HiTAC  of  Katsuki 
and  Hasegawa  [3  .  In  other  words,  FLOX  is  only  a  subset  of  HiTAC, 
and  hence  the  two  are  not  identical.  Unfortunately,  the  border  of 
FLOX  cannot  be  determined  from  refs.  1,  3,  and  4. 
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Cavaliere  and  de  Joannon  [1  analyzed  the  FLOX  and  HiTAC  and 
concluded  that  those  technologies  or  the  like  all  use  moderate  or 
intense  low-oxygen  dilution,  which  can  be  abridged  as  ‘MILD’. 
These  authors  defined  the  ‘MILD  combustion’  in  a  rigorous  and 
quantitative  fashion.  That  is,  the  occurrence  of  MILD  combustion 
requires:  (1 )  the  inlet  temperature  of  the  reactant  mixture  Tm  >  Ts j, 
and  (2)  the  allowable  maximum  temperature  rise  due  to  combus¬ 
tion  A T  <  TS[.  Based  on  these  conditions,  they  drew  a  very  simple 
Tin-AT  diagram  that  unambiguously  shows  different  combustion 
regimes  of  burning  CH4/O2/N2  in  a  WSR,  see  Fig.  1(c),  with  a  1.0  s 
residence  time  and  at  atmospheric  pressure;  they  termed  these 
regimes  as  the  feedback  combustion  (FC),  high  temperature  com¬ 
bustion  (HTC)  and  moderate  or  intense  low-oxygen  dilution  (MILD) 
combustion.  From  Refs.  1,3  ,  both  MILD  and  HTC  regimes  are 
subsets  of  HiTAC  of  Katsuki  and  Hasegawa  [3]  while  FC  is  somehow 
identical  to  PAC. 

Following  Cavaliere  and  de  Joannon  [1  ,  the  classification  of  the 
above  similar  technologies,  i.e.,  FLOX,  HiTAC  and  MILD  combustion, 
has  been  extensively  investigated  using  the  counter-flow  configu¬ 
rations  [11-14].  In  those  studies,  different  combustion  regimes  of 
the  diffusion  flames  had  been  well  achieved  and  evaluated  under 
various  non-premixed  combinations  of  hot  and  diluted  fuel/ 
oxidant  flow,  e.g.,  hot  oxidant  and  diluted  fuel  (HODF)  [11,12],  hot 
fuel  and  diluted  fuel  (HFDF)  [13],  and  also  hot  oxidant  and  diluted 
oxidant  (HODO)  14].  Note  that  all  those  previous  studies  cited 
[1-14]  on  the  FLOX,  HiTAC  or  MILD  combustion  were  performed  for 
the  highly-diluted  low-oxygen  cases.  However,  we  have  found  that 
the  quantitative  mathematical  conditions  of  T[n  >  TS[  and  AT  <  TS[  for 
MILD  combustion,  which  were  derived  by  Cavaliere  and  de  Joannon 
[1],  do  not  require  reactants  to  be  highly  diluted.  Based  on  this 
finding,  our  recent  work  15]  has  been  conducted  to  investigate  the 
combustion  regime  for  fuel-jet  flames  in  hot  co-flow  (JHC)  with  the 
oxygen  concentration  X02  =  0%— 100%  and  temperature 

TC0f  =  200  K-2400  K.  We  found  that  the  JHC  flames  can  be  also 
classified  into  three  regimes:  traditional  combustion  (TC),  high 
temperature  combustion  (HTC),  and  flameless  combustion  (FLC). 
The  FLC  regime  can  be  further  divided  into  three  distinct  sub-zones, 
i.e.,  MILD,  MILD-like  and  Quasi-MILD.  Here,  the  MILD  and  MILD-like 
combustion  both  satisfy  the  MILD  conditions  (TCOf  >  Ts j  and  AT  <  Ts 0 
of  Cavaliere  and  de  Joannon  [1]  while  the  quasi-MILD  combustion 
does  not.  Moreover,  the  MILD-like  combustion  occur  in  the  non- 
diluted  and  enriched  cases  where  X02  >21%  for  the  sufficiently 
high  co-flow  temperature.  Apparently,  the  MILD-like  combustion 
has  broken  the  ‘prerequisite’  of  low-oxygen  for  the  occurrence  of 
MILD  combustion.  Nevertheless,  it  is  not  clear  that  whether  this 
recent  finding  [15  is  still  appropriate  to  describe  the  combustion 
regimes  in  a  well-stirred  reactor  (WSR). 


On  the  other  hand,  de  Joannon  et  al.  [16,17  conducted  the 
analysis  of  methane  oxidation  in  a  well-stirred  reactor  (WSR)  to 
evaluate  the  diluted  MILD  combustion.  They  found  that  the  overall 
reaction  mechanism  of  CH4  oxidation  varies  with  different  WSR 
working  temperature  (Twsr)  16].  It  was  also  found  that  the 
competition  between  the  exoergic  oxidation  and  the  endothermic 
recombination  plays  an  important  role  under  MILD  condition  17]. 
As  Twsr  is  increased  to  1200  K,  the  C(2)  compounds  become  prev¬ 
alent  and  the  endothermic  recombination  is  favored  over  the 
oxidative  one.  However,  the  previous  studies  [16,17]  were  carried 
out  only  for  the  diluted  case,  leaving  the  non-diluted  and  enriched 
cases  of  combustion  open  for  further  investigations. 

Moreover,  elementary  chemical  pathways  of  fuel  oxidation 
differ  greatly  under  various  reacting  temperatures  [18  .  The  varia¬ 
tion  of  the  combustion  regime  may  result  in  different  reacting 
temperature  and  then  causes  the  pathway  change  of  the  fuel 
oxidation.  From  this  point,  Plessing  et  al.  [8]  carried  out  the 
instantaneous  measurements  of  temperature  and  OH  concentra¬ 
tion  in  a  recuperative  furnace.  They  found  that,  under  MILD  con¬ 
dition,  OH  is  homogeneously  distributed  to  a  relatively  larger 
region  and  its  concentration  in  the  combustion  zone  is  lower  than 
that  under  traditional  combustion  condition.  Compared  with  the 
traditional  combustion,  the  intermediate  products  such  as  OH  and 
CO  were  found  to  be  extremely  suppressed  when  operating  at  MILD 
condition,  which  has  also  reported  by  Dally  et  al.  19].  Using  a  fuel 
jet  in  hot  low-oxygen  coflow  (JHC)  device,  Dally  et  al.  19]  emulated 
the  MILD  combustion  and  deduced  that  there  would  be  difference 
of  the  chemical  pathways  between  the  highly-diluted  and  non- 
diluted  JHC  flames.  To  date,  however,  although  abundant  previous 
studies  cited  1-17]  have  been  performed  on  the  MILD  combustion, 
little  attention  has  been  paid  on  the  comparison  of  reaction  path¬ 
ways  of  fuel  oxidation  under  different  combustion  regimes. 

Accordingly,  to  fill  the  gaps  noted  above,  the  present  study  is 
designated  to  characterize  the  WSR  combustion  of  diluted,  non- 
diluted  and  enriched  CH4/02  mixture.  Based  on  the  ignition- 
extinction  S-curve,  the  extinction  and  self-ignition  temperatures 
(Tex  and  Ts 0  of  the  CH4/O2/N2  mixture  at  different  compositions  are 
systematically  investigated.  The  main  objectives  of  this  study  are 
twofold: 

(i)  To  quantitatively  classify  the  combustion  regimes  in  a  WSR 
over  the  full  range  of  oxygen  (0%-100%)  that  extends  the 
previous  understanding  of  MILD  combustion  derived  from 
Cavaliere  and  de  Joannon  1];  and 

(ii)  To  examine  the  overall  and  elementary  chemical  pathways  of 
CH4  oxidation  occurring  under  different  WSR  combustion 
regimes. 


Fig.  1.  Different  combustion  regimes,  (a)  Stability  limits  of  firing  CH4/air,  based  on  the  in-furnace  or  internal  parameters  [2  ,  where  FLOX  =  Flameless  Oxidation,  TC  =  Traditional 
Combustion,  and  UC  =  Unstable  Combustion;  (b)  auto-ignition  limit  (AIL)  and  forced-ignition  limit  (FIL)  of  C3H8  in  preheated  air  or  diluted  air  with  N2,  based  on  the  external 
parameters  [3,4];  (c)  Tjn  -  AT  locus  of  different  WSR  combustion  regimes  for  a  mixture  of  CH4/02/N2  at  atmospheric  pressure  and  a  residence  time  of  1  s,  where  FC  =  Feedback 
Combustion,  HTC  =  High  Temperature  Combustion  and  MILD  =  Moderate  or  Intense  Low-oxygen  Dilution  [1  . 
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Fig.  2.  Schematic  diagram  of  the  well  stirred  reactor  (WSR). 


2.  Computational  approach 

Present  numerical  simulations  of  the  premixed  CH4/02/N2 
combustion  in  a  WSR  are  performed  using  CHEMKIN  4.1  [20] 
software.  The  reactor  is  an  ideal  model,  where  reactants  can  be 
mixed  perfectly  before  reactions  take  place.  This  model  has  been 
used  to  study  many  aspects  of  combustion  such  as  flame  stabili¬ 
zation  and  NO*  formation.  The  reactions  in  the  WSR  take  place  in  a 
homogeneous  and  adiabatic  environment  and  may  be  approxi¬ 
mated  to  those  locally  occurring  at  a  spatial  point  in  practical  fur¬ 
naces.  Therefore,  spatial  distributions  of  species  usually 
encountered  in  industrial  combustors  are  not  considered  in  WSR. 
Fig.  2  shows  a  schematic  diagram  of  the  present  WSR  model.  On  the 
diagram,  V ,  P  and  r  are  the  volume,  pressure  and  mixture  residence 
time  in  the  reactor;  X,-  and  TWsr  are  the  mole  fraction  of  species  i 
and  the  mixture  temperature  in  the  reactor;  X*  and  T[n  denote  the 
mole  fraction  of  inlet  species  i  and  the  inlet  mixture  temperature; 
X°  and  Tout  indicate  the  mole  fraction  of  outlet  species  i  and  the 
outlet  mixture  temperature.  In  this  study,  V  and  P  are  taken  to  be 
67.4  cm3  and  1.0  atm.  The  residence  time  r  can  be  used  as  a  char¬ 
acteristic  parameter  of  the  reactor.  The  system  is  assumed  adiabatic 
and  so  the  heat  exchange  between  the  reactor  and  the  ambient  is 
set  zero.  The  total  mass  flow  rate  (M[n)  into  the  reactor  is  equal  to 
that  (Mout)  issuing  out  from  the  reactor  without  any  mass  exchange 
with  the  surrounding.  The  mass  flow  rate  M[n  =  Mout  is  related  to  V 
and  r  by  M[n  =  pVjr ,  where  p  is  the  density  of  the  gas  mixture  inside 


the  reactor.  To  summarize,  the  perfectly  premixed  CH4/02/N2 
mixture  at  T\n  flows  into  the  reactor,  CH4  reacts  with  02  when  T[n  is 
sufficiently  high  and  then  the  combustion  products  flow  out  from 
the  reactor.  The  mixture  residence  time  r  can  be  varied  by  changing 
the  mixture  mass  flow  rate  Mm  or  the  reactor  volume  V. 

All  calculations  use  the  detailed  mechanism  of  GRI-Mech  3.0 
[21  which  consists  of  53  species  for  a  total  of  325  reversible  re¬ 
actions.  GRI-Mech  3.0  is  an  optimized  mechanism  that  was 
designed  to  model  natural  gas  (mainly  CH4)  combustion  under  the 
conditions  of  temperature  varying  from  1000  K  to  2500  K,  pressure 
from  10  Torr  to  10  atm,  and  equivalence  ratio  from  0.1  to  5.0.  Note 
also  that  the  thermal  properties  involved  in  GRI-Mech  3.0  are 
ranged  from  200  K  to  3000  K.  The  previous  studies  have  shown  its 
ability  in  modeling  air-diluted  [22],  CH4/air  [23],  and  oxy-enriched 
diffusion  flames  [24  . 

3.  Modeling  validation 

To  validate  the  present  modeling,  we  model  the  CH4/air  com¬ 
bustion  in  the  WSR  under  the  same  conditions  as  those  used  in  the 
experiments  of  Steven  and  Joseph  25] .  The  temperature  of  the  inlet 
CH4  and  air  is  298  K.  The  CH4/air  mixture  is  burned  over  the 
equivalence  ratio  from  0.53  to  1.63  and  the  average  residence  time 
of  the  mixture  in  the  reactor  is  0.0085  s. 

Fig.  3  compares  the  present  calculations  with  the  experimental 
data  and  the  previous  predictions  of  Steven  and  Joseph  [25]. 
Obviously,  the  present  WSR  and  previous  numerical  results  both 
predict  well  Twsr  and  X02  over  the  broad  equivalence  ratio  range. 
Consistent  with  the  measurements,  the  predictions  show  a  slight 
rise  of  X02  when  cp  >  1.4.  Apparently,  at  the  fuel  rich  conditions  of 
(p  >  1.4,  Twsr  is  relatively  low  (see  Fig.  3(a))  that  the  oxidation  of  CH4 
cannot  process  completely  over  the  residence  time.  Moreover, 
although  the  predictions  of  Xco2  are  less  satisfactory,  their  variable 
trend  is  similar  to  the  measurement.  Accordingly,  the  present  WSR 
model  with  GRI-Mech  3.0  is  effective  in  simulating  the  WSR  com¬ 
bustion  of  CH4/air  and  hence  can  be  used  in  the  subsequent 
sections. 


(a) 

0.40  0.75  1.10  1.45  1.80 


Fig.  3.  Comparisons  between  the  predictions  and  the  measurements  [25]  of  (a)  Twsr.  (b)  X02  and  (c)  XCo2 
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Fig.  4.  WSR  working  temperatures  (Twsr)  versus  (a)  the  inlet  mixture  temperature  (Tjn)  and  (b)  the  residence  time  (t)  for  the  CH4/O2/N2  mixture  at  6.7%  O2,  3.3%  CH4  and  90%  N2. 


4.  Ignition  and  extinction  of  WSR  combustion 

Fig.  4  plots  Twsr  versus  T[n  with  T[n  =  300  K-3000  I<  and  TWsr 
versus  r  with  the  residence  time  (10  5  s  ~  10  s)  for  the  reactant 
mixture  (6.7%  O2  +  3.3%  CH4  +  90%  N2).  In  Fig.  4,  the  well-known 
ignition-extinction  S-curves  26]  are  illustrated.  The  upper  and 
lower  branches  of  the  S-curve  represent  the  states  of  strong 
burning  and  weak  oxidation  whereas  the  middle  branch  of  the  S- 
curve  (dashed  line)  is  the  unstable  state  that  cannot  be  obtained 
experimentally.  In  Fig.  4(a),  T\n  =  Tex  at  the  upper  turning  point  is 
the  maximum  temperature  of  extinction  of  the  strong  burning 
flame,  i.e.,  ‘extinction  temperature’,  while  T\n  =  Ts j  at  the  lower 
turning  point  is  considered  as  the  minimum  temperature  of  reac¬ 
tant  self-ignition,  often  termed  as  ‘self-ignition  temperature’. 
Similarly,  in  Fig.  4(b),  the  turning  point  (r  =  rex)  of  the  upper  branch 
is  the  maximum  residence  time  below  which  the  flame  extinction 
occurs  whereas  the  lower  turning  point  of  r  =  rSi  is  considered  as 
the  minimum  residence  time  of  the  self-ignition. 

In  Fig.  4(a),  as  T\n  is  increased  to  Ts j,  the  mixture  reaction  occurs 
spontaneously  and  Tws r  jumps  from  the  lower  branch  up  to  the 
upper  branch;  then  Twsr  increases  along  the  upper  branch  as  T\n  is 
further  increased.  Fig.  4(a)  also  demonstrates  that  both  Ts j  and  ?ex 
climb  as  r  is  shortened.  This  can  be  easily  understood  here.  A 


Fig.  5.  WSR  working  temperatures  (Twsr)  versus  the  inlet  mixture  temperature  (Tin) 
with  r  =  1.0  s  for  three  different  stoichiometric  mixture  compositions  of  CH4/O2/N2, 
i.e.,  6.7%  02  +  3.3%  CH4  +  90%  N2, 19%  02  +  9.5%  CH4  +  71.5%  N2,  66.7%  02  +  33.3%  CH4 
and  0.67%  02  +  0.33%  CH4  +  99%  N2. 


shorter  r  causes  less  heat  released  in  the  reactor  and  hence  requires 
a  higher  Tm  for  the  reaction  maintaining.  Fig.  4(b)  shows  that,  at 
Tin  =  600  K,  the  WSR  mixture  may  be  not  ignited  by  increasing  the 
residence  time.  In  contrast,  when  Tm  is  escalated,  the  residence 
time  for  the  ignition  is  reduced  significantly.  This  should  result 
from  the  fact  that  the  increase  of  Tm  makes  the  reactants  more 
active,  which  hence  allows  the  rex  and  rSi  to  decrease.  In  addition,  as 
Tin  is  increased,  the  difference  between  Ts \  and  Tex  narrows  and 
finally  disappears. 

Moreover,  Cavaliere  and  de  Joannon  [1  analyzed  the  different 
combustion  regimes  of  CH4/O2/N2  combustion  in  a  WSR  with  the 
residence  time  r  =  1.0  s  and  at  atmospheric  pressure  (1  atm).  For 
consistency  with  their  investigation,  this  study  also  uses  this  value 
of  r  for  calculations  presented  hereafter. 

Fig.  5  illustrates  the  Tw sr  against  T\n  with  r  =  1  s  for  four  mixture 
compositions,  i.e.,  66.7%  O2  +  33.3%  CH4  (O2  +  CH4),  19%  O2  +  9.5% 
CH4  +  71.5%  N2  (air  +  CH4),  6.7%  02  +  3.3%  CH4  +  90%  N2  (diluted 
oxidant  +  CH4),  and  0.67%  O2  +  0.33%  CH4  +  99%  N2  (extremely 
diluted  oxidant  +  CH4).  Results  demonstrate  that,  as  X^2  is 
increased,  the  middle  branch  shifts  up  to  higher  values  of  Tmi  which 
was  also  reported  in  Refs.  [1,16].  That  is,  both  TS[  and  Tex  increase 
when  adding  more  N2  into  the  reactant  mixture.  Moreover,  similar 
to  its  dependence  on  T\n  and  r  (Fig.  4),  the  difference  between  Ts j 
and  Tex  narrows  with  increasing  XjJ,  and  nearly  vanishes  at 


Fig.  6.  WSR  working  temperatures  (Twsr)  at  the  upper  steady  state  as  a  function  of  the 
inlet  temperature  (Tjn)  for  stoichiometric  CH4/O2/N2  mixture. 
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Fig.  7.  WSR  temperature  increment  in  the  upper  steady  state  (AT  =  Tws r  -  Tin)  versus 
the  inlet  temperature  (Tjn)  for  the  stoichiometric  CH4/02/N2  mixture.  On  the  map, 
NC  =  no  combustion;  TC  =  traditional  combustion;  HTC  =  high  temperature  com¬ 
bustion;  MILD  =  moderate  or  intense  low-oxygen  dilution. 

XjJ,  =  99%  or  higher.  This  is  also  consistent  with  the  previous  ob¬ 
servations  for  the  low-oxygen  combustion  under  the  WSR  [27  or 
the  counter-flow  diffusion  conditions  [28  . 

5.  Classification  of  WSR  combustion 

Fig.  6  plots  Twsr  ^  the  upper  steady  state  against  Tm  for  the 
stoichiometric  CH4/O2  mixture  with  XjJj  ranging  from  0%  to  99%. 
Obviously,  TWsr  decreases  as  Xj^  is  increased  at  a  fixed  T[n.  This  is 
because,  at  a  higher  XjJj ,  the  amount  of  the  inlet  reactants  (CH4/O2) 
is  less  and  thus  less  heat  per  unit  the  mixture  is  released,  which 
consequently  leads  to  a  lower  Twsr-  For  each  mixture  composition, 
Twsr  increase  monotonically  with  T\n.  The  reason  is  that  the  re¬ 
actants  initially  carry  more  enthalpy  at  a  higher  T[n.  Similar  obser¬ 
vations  were  also  reported  in  refs.  1  and  16,  but  their  discussions 
were  all  limited  for  X^2  >  71.5%.  The  Tin-TWSR  map  can  be  easily 
transformed  to  a  Tin-AT one,  where  AT  is  the  working  temperature 
increment,  i.e.,  AT  =  TWsr  -  T\n.  Since  TS[  varies  with  the  CH4/O2/N2 


(a)  x*N  =  0% 

'  v2 


Table  1 

The  classification  of  different  regimes  for  CH4/N2/02  combustion  in  a  WSR. 


Combustion  regime 

Inlet  conditions 

Working  conditions 

Traditional  combustion  (TC) 

Tn  <  Tsi 

AT  >  Tsi 

High  temperature  combustion  (HTC) 

Tin  >  Tsj 

AT  >  Tsi 

Flameless  combustion  (FLC) 

Quasi-MILD  combustion 

Tin  <  Tsi 

AT  <  Tsi 

MILD  +  MILD-like  combustion 

Tin  >  Tsj 

AT  <  Tsi 

composition,  the  abscissa  and  ordinate  of  the  Tin-AT  plot  are 
normalized  using  Tsl. 

As  shown  in  Fig.  7,  for  a  fixed  T[ni  AT  decreases  rapidly  with 
increasing  XjJ,  due  to  the  reduction  of  Tw sr-  More  importantly,  for 
any  mixture  composition,  AT  decreases  with  T\n,  which  can  be 
explained  here.  The  increase  of  Tm  leads  to  a  higher  WSR  working 
temperature,  which  is  more  beneficial  to  the  dissociation  of  CO2 
and  H2O  18].  Evidently,  a  check  of  the  outlet  product  composition 
reveals  that  an  increase  in  Tm  yields  more  unburnt  combustible 
gases  such  as  CO  and  H2,  see  Fig.  8.  In  other  words,  as  T[n  is 
increased,  less  complete  combustion  takes  place  in  the  reactor  and 
thus  less  reaction  heat  is  released  over  the  same  residence  time, 
hence  reducing  AT.  In  addition,  the  heat  capacities  of  flue  gases 
(e.g.,  CO2,  O2,  N2,  and  H2O)  are  expected  to  increase  with  temper¬ 
ature  18],  and  so  for  higher  Fin>  the  same  heat  enables  the  flue 
gases  to  have  a  lower  temperature  rise.  To  summerize,  both  effects 
of  the  above  chemical  and  physical  factors  due  to  increasing  Tm 
result  in  the  reduction  of  AT. 

Following  the  method  used  by  Cavaliere  and  de  Joannon  1]  for 
the  WSR  combustion,  a  vertical  line  at  Tin/TSi  =  1  and  a  horizontal 
line  at  AT/TSi  =  1  are  drawn  to  divide  the  map  into  several  regions, 
see  Fig.  7.  That  is,  the  WSR  combustion  can  be  classified  as  the 
following  four  regimes,  i.e.,  no  combustion  (NC),  traditional  com¬ 
bustion  (TC),  high  temperature  combustion  (HTC),  and  flameless 
combustion  (FLC),  where  FLC  comprises  the  MILD,  or  MILD-like, 
and  Quasi-MILD  sub-regimes.  The  criterion  for  the  various  com¬ 
bustion  regimes  is  summarized  in  Table  1. 

In  the  NC  region,  T\n  is  less  than  Tex  and  the  reactants  are  so 
diluted  that  the  released  heat  by  the  reaction  cannot  sustain  the 
oxidation  process.  Hence,  the  combustion  is  extinct.  The  TC  regime 
is  located  in  the  upper  left  of  the  map,  where  T\n  is  low  and  the 
working  temperature  increment  is  high,  reflecting  the  typical 


(b)  X*  =  89.8% 

/V2 


Fig.  8.  Outlet  product  concentrations  XHz ,  Xco2  and  Xco  versus  the  inlet  temperature  ( 7% )  for  two  stoichiometric  CH4/02/N2  mixture  compositions:  (a)  66.7%  02  +  33.3%  CH4  and  (b) 
6.7%  02  -1-  3.3%  CH4  +  90%  N2.  On  the  map,  X£o  is  the  theoretical  C02  concentration  of  the  flue  gases;  NC  =  no  combustion;  TC  =  traditional  combustion;  HTC  =  high  temperature 
combustion;  MILD  =  moderate  or  intense  low-oxygen  dilution. 
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Fig.  9.  Classification  of  the  premixed  stoichiometric  CH4/O2/N2  combustion  in  a  WSR. 
On  the  map,  Ts j  and  Tex  are  the  self-ignition  and  extinction  temperatures  while 
AT  =  TWsr  -  Tin.  NC  =  no  combustion;  TC  =  traditional  combustion;  HTC  =  high 
temperature  combustion;  MILD  =  moderate  or  intense  low-oxygen  dilution. 

normal  flames.  As  T[n  is  increased,  the  combustion  state  is  trans¬ 
formed  to  the  HTC  regime  where  the  combustion  features  differ 
from  those  of  the  TC.  The  lower  region  of  the  map  belongs  to  the  FLC 
because  the  temperature  rise  is  sufficiently  small  at  AT  <  Tsl  so  that 
the  flame  might  be  invisible.  More  specifically,  on  the  right  of  the 
FLC  regime  is  the  MILD  or  MILD-like  combustion,  where  the  high  T[n 
of  the  reactants  promotes  and  sustains  the  reaction  process  with  a 
temperature  rise  smaller  than  Ts j,  i.e.,  T\n  >  Ts j  and  A T  <  Ts j.  On  the 
left,  a  triangular  area  represents  a  transitional  region  from  the  NC  to 
the  TC,  HTC  and  MILD  or  MILD-like  combustion.  This  transitional 
zone  is  featured  by  Tm  <  Tsl  and  A T  <  Ts j.  In  this  zone,  although  Tm  is 
not  high  enough  to  spontaneously  ignite  the  reactants  or  cause  fast 
combustion  reactions,  the  accumulation  of  the  input  heat  energy 
carried  by  the  preheated  initial  mixture  and  the  heat  released  from 
exothermic  reactions  can  sustain  slow  or  lmild’  combustion.  Based 
on  these  characteristics,  this  transitional  zone  is  named  as  ‘Quasi- 
MILD’  in  the  FLC  regime.  Note  that  similar  observations  to  the 
transitional  zone  were  also  reported  in  Refs.  12,13,15]. 

Fig.  8  shows  the  concentrations  of  the  outlet  unburnt  combus¬ 
tible  gases  (CO  and  H2)  versus  the  inlet  temperature  Tjn.  Quite 
obviously,  as  T\n  is  increased,  the  emissions  of  CO  and  H2  grow 


steadily.  That  is,  at  a  higher  T[n,  the  conversion  from  CH4  to  CO2 
becomes  less  complete.  Fig.  8(a)  shows  that  the  CO2  emission  (XCo2 ) 
is  far  less  than  the  theoretical  one  (X^  )  for  the  non-dilution  case. 
This  indicates  that  the  conversion  rate  from  CH4  to  CO2  is  very  small 
no  matter  at  which  regime  the  combustion  occurs.  In  contrast,  for 
the  diluted  case  (Fig.  8(b)),  the  initial  CH4  is  nearly  completely 
oxidized  into  CO2,  i.e.,  XCo2  ~^co2  =  ^ch4  when  Tm  <  1300  K,  which 
means  the  high  conversion  rate  from  CH4  to  CO2  under  quasi-MILD 
or  MILD  combustion  regimes.  However,  when  T\n  >  1300  K,  XC02 
begins  to  decrease  with  further  increasing  Tjn.  Hence,  the  conversion 
from  CH4  to  CO2  becomes  less  complete  at  high  inlet  temperatures 
even  under  MILD  or  MILD-like  conditions.  From  this,  although  MILD 
and  MILD-like  combustion  can  well  satisfy  the  mathematical  defi¬ 
nition  1  by  increasing  T\n  for  any  reactant  composition,  the 
occurrence  of  these  regimes,  especially  for  non-diluted  and  enriched 
cases,  is  achieved  at  the  expense  of  reduction  in  conversion  rate  from 
CH4  to  CO2.  Therefore,  when  referred  to  the  mathematical  definition 
[  1  of  MILD  combustion,  the  conversion  rate  from  CH4  to  C02  should 
also  be  considered  in  the  indentification  of  MILD  combustion.  Be¬ 
sides,  considering  the  adiabatic  conditions  of  the  present  WSR  cal¬ 
culations,  the  reaction  heat  released  cannot  be  transferred  outside 
via  the  reactor  wall.  However,  it  is  important  to  note  that  this  situ¬ 
ation  will  change  if  the  WSR  is  replaced  by  a  practical  combustor 
such  as  a  heating  furnace.  In  the  latter,  the  reaction  heat  is  contin¬ 
uously  transferred  away  and  hence  the  reacting  temperature  drops 
along  the  reaction  process,  finally  terminating  the  dissociation  of 
CO2  and  H2O.  Therefore,  the  fuel  will  surely  be  burnt  out  for  (f>  <  1 
with  sufficiently  long  reacting  or  residence  time. 

With  the  aid  of  Tin-AT  map,  various  combustion  regimes  can  be 
identified  under  different  combustor  configurations,  as  shown  in 
the  numerous  existing  publications  [1,11—14  .  However,  the  short¬ 
coming  of  using  Tin-AT  map  is  obvious  that  one  cannot  determine 
the  combustion  regime  a  priori  because  the  operative  conditions 
are  given  by  Tjn,  reactant  concentrations  (X/),  and  q>,  rather  than  AT. 
Consequently,  the  Tin-AT  map  is  not  available  to  directly  determine 
which  regime  the  combustion  will  be  evolved  under  given  condi¬ 
tions  (Tin,  X/,  and  <p).  To  address  this  issue,  we  transform  the  Tjn— AT 
map  of  Fig.  7  into  a  T[n—X^  map,  shown  in  Fig.  9.  On  the  map,  the 
four  aforementioned  combustion  regimes  are  clearly  indentified 
using  the  three  lines,  i.e.,  Tin  =  Tex,  Tin  =  TSi  and  AT  =  TSi.  Fig.  9 
demonstrates  that  the  mathematical  definition  of  MILD  combus¬ 
tion  [1]  can  be  well  satisfied  as  long  as  Tin  is  sufficiently  high.  In 
other  words,  the  MILD  and  MILD-like  combustion  as  a  whole  can 
occur  for  all  cases  even  without  reactant  dilution.  For  instance, 
MILD  or  MILD-like  combustion  can  occur  in  the  cases  with 


(a)^  —  0.5  (b)  (/) -  1  (c)  (j)  2 


Fig.  10.  Classifications  of  the  premixed  CH4/O2/N2  WSR  combustion  at  (a)  <p  =  0.5,  (b)  cp  =  1.0,  and  (c)  cp  =  2.0.  On  the  map,  NC  =  no  combustion;  TC  =  traditional  combustion; 
HTC  =  high  temperature  combustion;  MILD  =  moderate  or  intense  low-oxygen  dilution. 
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Fig.  11.  Schematic  3D  diagram  showing  the  combustion  regimes  expressed  in  terms  of 
(p,  Tin  and  X^  .  On  the  diagram,  NC  =  no  combustion;  TC  =  traditional  combustion; 
HTC  =  high  temperature  combustion;  MILD  =  moderate  or  intense  low-oxygen 
dilution. 

Tin  >  2520  K  and  X^2  =  0%  and  those  with  T[n  >  945  K  and  X^2  =  90%. 
Moreover,  through  Fig.  9,  one  can  directly  determine  the  final 
combustion  regime  from  any  given  operative  conditions  (Tin,  X/) 
while  it  is  hard,  at  least  inconvenient,  to  make  such  estimations 
from  the  Tin-AT  map.  On  the  other  hand,  the  overall  combustion 
regime  classification  depends  little  on  the  combustor  configura¬ 
tions,  because  the  same  regimes  can  be  also  classified  in  the  pre¬ 
vious  diffusion  systems  [12,13,15]. 


In  fact,  the  local  equivalence  ratio  inside  a  furnace  cannot  be 
always  unity.  Accordingly,  to  emulate  the  combustion  inside  a 
practical  furnace,  the  effect  of  <p  is  also  examined.  Fig.  10  indicates 
the  classifications  of  different  combustion  regimes  at  <p  =  0.5  (fuel 
lean),  1.0  (stoichiometric)  and  2.0  (fuel  rich).  As  expected,  for  the 
three  values  of  <p,  the  three  maps  are  overall  similar,  but  the  loca¬ 
tions  and  sizes  of  the  regimes  change.  Apparently,  the  line  of 
Tin  =  rSi  is  lower,  whereas  that  of  A T  =  Ts j  is  higher,  for  the  stoi¬ 
chiometric  mixture  than  for  the  other  two  cases.  This  is  because  the 
fuel  can  be  more  easily  ignited  and  burned  out  under  the  stoi¬ 
chiometric  condition.  Also,  under  the  same  dilution  level,  Tw sr  for 
<p  =  1.0  should  be  higher  than  those  of  <p  =  0.5  and  <p  =  2.0,  so  that  a 
higher  Tm  is  required  to  satisfy  the  MILD  combustion  condition 
AT  =  TwSR-Tin  <  Tsi. 

Based  on  the  results  of  Fig.  10  and  those  for  other  values  of  <p,  the 
classification  for  the  CH4/O2/N2  combustion  can  be  illustrated 
schematically  in  a  three-dimensional  plot  (Tin  -  <p  -  XjJj ),  see  Fig.  11. 
It  is  demonstrated  that  the  extinction  limit  is  extended  to  higher 
xn2  as  Tin  is  increased  and/or  cp  *  1.  The  Tjn  —  Tgx  plane  represents 
the  critical  operative  conditions  beyond  which  no  combustion  oc¬ 
curs.  The  reaction  region  is  then  divided  into  two  regions  by  the 
Tin  =  TSi  plane,  below  which  the  reactant  self-ignition  cannot  take 
place.  The  addition  plane  of  the  AT  =  Tsj  enables  the  WSR  com¬ 
bustion  to  be  subdivided  into  five  regimes,  i.e.,  no  combustion  (NC), 
traditional  combustion  (TC),  high  temperature  combustion  (HTC), 
Quasi-MILD  combustion,  MILD  or  MILD-like  combustion.  Here,  the 
MILD  or  MILD-like  and  Quasi-MILD  combustion  constitute  the 
flameless  combustion  (FLC).  From  Fig.  11,  one  can  easily  identify  the 
combustion  regime  for  any  given  conditions  (temperature,  re¬ 
actants  equivalence  ratio,  and  dilution  ratio).  Moreover,  Medwell 
and  Dally  29]  found  that  the  fuel  type  does  not  have  a  significant 
effect  on  the  combustion  reaction  characteristics.  Therefore,  when 
referred  to  other  fuels,  the  combustion  regimes  shown  in  Fig.  11 
should  also  apply. 
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Fig.  12.  Conversions  of  the  initial  C  and  H  to  CO,  C02,  C2H2,  H20  and  H2  as  functions  of  7WSr  at  r  =  1.0  s  and  <p  =  4.0.  Solid  black  lines  represent  the  C  and  H  conversions  whereas  the 
marked  areas  denote  the  contributions  of  C  and  H  conversions  to  the  productions  of  CO,  C02,  C2,  H20  and  H2  at  three  different  N2  dilution  levels. 
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Fig.  13.  Conversions  of  the  initial  C  and  H  to  CO,  CO2,  H2,  H2O,  H  and  OH  as  functions  of  TWsr  at  cp  =  1.0  and  r  =  1.0  s.  The  three  different  N2  dilution  levels  are  indicated  on  the  plots. 


6.  Characteristics  of  WSR  combustion 

6.2.  Overall  mechanism  of  CH 4  oxidation 

Following  the  above  discussion,  it  is  aware  that  the  transition 
between  the  combustion  regimes  and  different  Tw sr  can  be  ach¬ 
ieved  by  varying  the  inlet  reactant  conditions,  i.e.,  <p,  T[n  and  Xj^. 
One  would  then  ask  whether  the  products  composition  remain  the 


same  under  different  combustion  regimes.  It  is  hence  worthy  to 
investigate  the  effects  of  c Pi  T'm  and  X^2  on  the  overall  reaction 
mechanisms  of  CH4  combustion  when  varying  combustion  regimes. 

Fig.  12  presents  the  contributions  of  the  initial  carbon  (C)  and 
hydrogen  (FI)  atom  to  the  productions  of  CO,  CO2,  C2H2,  H2O  and  H2 
as  functions  of  Tw sr  at  r  =  1.0  s  and  cp  =  4.0.  As  claimed  in  Ref.  16], 
C2H2  is  the  most  important  one  of  the  C2  compounds  (C2H6,  C2H5, 
C2H4,  C2H3,  and  C2H2)  formed  through  recombination.  Thus,  only 
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Fig.  15.  Products  of  CH4  oxidation  dependence  on  cp  and  TWSr. 


C2H2  are  presented  here.  The  solid  lines  on  the  plots  represent  the 
ratios  between  the  converted  and  inlet  moles  of  atomic  C  (left)  and 
H  (right).  The  colored  areas  denote  the  C  and  H  yields  of  the  main 
produced  species  (C  CO,  CO2,  C2H2  and  H  H2O,  H2,  C2H2).  In 
addition,  the  curves  of  Fig.  12(a)  and  (b)  start  at  TWsr  =  1242  K, 
where  Tex  =  1050  I<  is  located  for  cp  =  4.0  andXj^  =  92%.  In  contrast, 
the  starting  points  of  the  curves  in  Fig.  12(c-f)  are  TWSr  =  1052  K 
and  1520  K  at  Tin  =  300  K  for  X^2  =  56%  and  0%. 

For  Twsr  <  1300  K,  Fig.  12(a-d),  CO,  H2,  and  H2O  are  the  main 
products  with  a  small  amount  of  CO2  and  almost  no  C2H2.  For 
1300K  <  Twsr  <  2000  K,  Fig.  12  demonstrates  that  the  concentra¬ 
tions  of  CO  and  H2  increase  rapidly  while  those  of  CO2  and  H2O 
decrease  and  C2H2  becomes  important.  This  is  because  the 
recombination  channel,  through  which  CH4  can  be  partly  recom¬ 
bined  to  the  C(2)  compound,  becomes  more  important  than  the 
oxidation  channel  in  this  range  of  TWsr  [16].  Then,  a  further  in¬ 
crease  of  Twsr  leads  all  the  C  and  H  to  be  conversed  to  CO  and  H2, 
and  hence  CO2,  H2O  and  C2H2  further  decrease  to  zero.  The  present 
results,  Fig.  12(a),  are  nearly  the  same  as  reported  in  Ref.  [16].  The 
authors  stated  that,  for  the  CH4/O2/N2  mixture  at  <p  =  4  and 
XjJ,  =  85%,  the  CH4  reaction  mechanisms  vary  greatly  with  TWsr- 
Based  on  Fig.  12,  the  most  representative  overall  mechanisms  for 
the  CH4  oxidation  can  be  expressed  as 

(i)  for  Twsr  <  1300  K:  CH4  +  02  -►  CO  +  H2  +  H20; 

(ii)  for  1300  <  TWsr  <  2000  K:  CH4  +  02  —  CO  +  H20  +  H2  +  C2H2; 

(iii)  for  Twsr  >  2000  K:  CH4  +  02  —  CO  +  H2. 

Figs.  13  and  14  show  the  ratios  of  C  and  H  yielding  the  products 
(i.e.,  CO,  CO2,  H2,  H2O,  H  and  OH)  for  <p  =  1.0  and  0.5.  In  these  figures, 
the  solid  lines  denote  the  ratios  between  the  converted  and  inlet 
moles  of  atomic  C  (left)  and  H  (right).  ForX^  =  92%,  the  conversion 
curves  start  at  TWsr  ~  500 1<  and  400 1<  for  cp  =  1.0  and  0.5.  For  other 
cases,  CH4  is  completely  consumed  and  hence  the  H  and  C  con¬ 
version  ratios  are  all  equal  to  1  (see  Figs.  13  (c-f)  and  14  (c-f)).  The 
colored  areas  represent  the  ratios  of  C,  H  atom  conversion  to  the 
main  products  (C  CO,  CO2  and  H  H2O,  H2,  H,  OH).  In  addition, 
very  little  C2  compounds  are  produced  for  these  cases,  and  hence 
they  are  not  given  on  the  plots.  Similar  as  that  for  <p  =  4.0,  the 
yielding  ratios  of  C  and  H  to  the  products  also  appear  to  be  a 
function  of  Twsr  for  the  two  values  of  cp  at  any  diluting  level.  CO2 
and  H20  are  found  to  be  the  main  products  for  TWsr  <  2000  K,  as 
shown  in  Figs.  13(a-d)  and  14(a-d).  For  TWsr  >  2000  K,  CO,  H2,  H, 


and  OH  become  important,  while  the  ratios  of  C  and  H  converted  to 
CO2  and  H2O  decrease.  In  these  cases,  the  global  mechanisms  of  the 
CH4  oxidation  can  be  proposed  as: 

(iv)  for  Twsr  <  2000  K:  CH4  +  02  —  C02  +  H20; 

( v)  for  Twsr  >  2000  K:  CH4  +  02  -►  CO  +  C02  +  H20  +  H2  +  H  +  OH. 

Following  the  overall  reactions  given  in  (i)  ~  (v),  Fig.  15  indicates 
the  impacts  of  TWsr  and  </>  on  the  products  of  CH4  oxidation.  For 
fuel-rich  conditions,  the  increase  of  Twsr  leads  to  three  composi¬ 
tions  of  the  products:  (i)  CO  +  H2  +  H2O  +  CO2,  (ii) 
CO  +  C2H2  +  H2  +  H20  and  (iii)  CO  +  H2.  In  contrast,  no  C(2) 
compound  can  be  produced  in  the  stoichiometric  or  fuel-lean 
conditions.  The  increase  of  TWsr  hence  only  results  in  the  two 
compositions  of  the  products:  (iv)  CO2  +  H2O  and  (v) 
CO  +  H2  +  H2O  +  CO2  +  OH  +  H.  Moreover,  according  to  Fig.  11,  the 
combustion  regime  transition  can  be  achieved  by  changing  Tm. 
Thus,  the  WSR  combustions  for  the  cases  of  Figs.  12-14  have  all 
experienced  the  regime  change.  Of  interest,  the  fuel  conversion  and 
its  corresponding  contributions  to  the  products  appear  to  be  in¬ 
dependent  of  any  combustion  regime.  Hence,  the  combustion 
regime  has  little  influence  on  the  overall  CH4  reaction  paths.  In 
other  words,  the  global  kinetic  mechanism  for  CH4  oxidation 
should  depend  mainly  on  Tw sr  and  <p  rather  than  the  initial  dilution 
level  and  the  combustion  regime. 

6.2.  Elementary  reactions  of  oxidizing  CH4 

Based  on  the  WSR  calculations,  the  elementary  chemical  path¬ 
ways  of  the  conversion  from  CH4  to  CO2  for  cases  1-6  ( Table  2)  are 
displayed  in  Fig.  16.  On  the  diagrams,  arrows  represent  elementary 
reactions  with  the  primary  reactant  species  at  the  tail  and  the 
primary  product  species  at  the  head.  Other  reactant  species  are 
shown  aside  of  the  arrows  and  the  corresponding  reactions 
(marked  in  accordance  with  GRI-Mech  3.0)  are  also  indicated.  The 
arrow  width  gives  a  visual  indication  of  the  relative  importance  of  a 
particular  reaction,  where  the  numbers  in  parentheses  quantita¬ 
tively  show  the  reaction  rates.  For  instance,  2.4-7  means  2.4  x  10-7 
(gmol/cm3  s).  In  the  diagrams,  unimportant  reactions  are  ignored. 
In  other  words,  only  those  pathways  with  reaction  rates  higher 
than  1  x  10-8  gmol/cm3  s  and  1  x  10~9  gmol/cm3  s  for  cases  1-5 
and  case  6,  are  shown. 

The  main  pathways  of  CH4/O2  mixture  reactions  at  different 
operative  conditions  can  be  summarized  briefly  as  follows: 

i)  For  Cases  1-3,  6: 

As  shown  in  Fig.  16(a~c,  f),  CH4  converses  mainly  through  the 
central  path  to  CO2,  with  several  side  paths  originating  from  the 
methyl  (CH3)  radical.  For  the  central  path,  CH4  molecule  firstly 
undergoes  H  abstraction  due  to  the  attack  by  OH,  O,  and  H  radicals 
and  produces  the  methyl  radical  (CH3).  The  CH3  then  combines 
with  an  oxygen  atom  to  form  formaldehyde  (CH2O),  which,  in  turn, 
is  attacked  by  OH  and  H  radicals  to  produce  the  formyl  radical 


Table  2 

Operative  conditions  for  CH4  reaction  analysis  under  different  combustion  regimes. 


Case 

Mixture  composition 

Tm 

Twsr 

Combustion  regime 

1 

33%  CH4  +  67%  02 

700  K 

3090  K 

TC 

2 

33%  CH4  +  67%  02 

1500  K 

3189  K 

HTC 

3 

33%  CH4  +  67%  02 

2700  K 

3338  K 

MILD-like 

4 

92%  N2  +  2.7%  CH4  +  5.3%  02 

700  K 

1312  K 

Quasi-MILD 

5 

92%  N2  +  2.7%  CH4  +  5.3%  02 

1500  K 

2052  K 

MILD 

6 

92%  N2  +  2.7 %  CH4  +  5.3%  02 

2700  K 

2867  K 

MILD 

CH4 


CH4 


M 

R12  (1.7-8) 


^OH  R99  (1.5-7) 

CC>2 


Fig.  16.  Reaction  pathways  for  the  non-diluted  (a~c:  cases  1-3)  and  diluted  (d~f:  cases  4~6)  CH4  combustion  for  r  =  1.0  s  in  a  WSR,  under  different  combustion  regimes,  i.e.,  case  1:  Traditional  combustion  (TC),  case  2:  high  temperature 
combustion  (HTC),  case  3:  MILD-like  combustion,  case  4:  Quasi-MILD  combustion  and  case  5,  6:  MILD  or  MILD-like  combustion.  Note  that,  the  parenthetical  numbers  denote  the  reaction  rates,  e.g.,  2.4-7  meaning  2.4  x  10~7  (gmol / 
cm3  s). 
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(HCO).  Next,  the  HCO  reacts  with  H2O  and  O2  to  form  CO.  Finally,  CO 
is  converted  to  CO2,  mainly  reacting  with  OH. 

For  the  side  paths,  all  start  from  CH3.  On  the  left,  CH3  reacts  with 
OH  radical  to  form  methylene  radical  in  two  possible  electronic 
configurations  (CH2(s)  and  CH2).  CH2(s)  reacts  with  H2O  and  N2  and 
is  translated  to  CH2,  which  forms  CH20,  HCO  and  then  CO;  note  that 
CH2(s)  only  reacts  with  H2O  in  non-diluted  cases  because  the  initial 
reactant  mixture  does  not  contain  any  N2.  Or  CH2  is  dehydro¬ 
genated  to  form  CH  and  C,  which  is  then  transformed  to  HCO  and 
CO.  Moreover,  there  are  two  different  paths  on  the  right  side:  i.e., 
CH3  is  either  converted  to  CH2OH  and  then  forms  CH2O,  HCO  and 
CO  or  is  directly  transformed  to  CO  and  then  C02. 

ii)  For  Case  5: 

It  differs  only  slightly  from  cases  1-3  and  6.  The  two- 
dehydrogenation  reactions  CH2  -►  CH  — ►  C  — ►  CO  and 
CH2  -►  CH  ->  HCO  CO  [green  arrows  in  Fig.  16  (a~c,  f)]  are  found 
to  disappear. 

iii)  For  Case  4: 

Apart  from  all  pathways  shown  in  Fig.  16(a~c  and  e~f),  those 
paths  that  are  unimportant  or  even  absent  for  cases  1-3,  5  and  6 
become  prominent  for  case  4.  As  demonstrated  in  Fig.  16(d),  red 
arrows  show  the  new  important  pathways  at  relative  low  reacting 
temperature  and  complement  the  diagram  for  the  high  tempera¬ 
ture  cases.  In  the  path  map,  several  interesting  pathways  appear. 
Firstly,  CH3  is  found  to  strongly  react  backward  to  CH4.  Then,  a  new 
route  from  CH3  to  CH2O  appears,  along  which  CH3  reacts  with  OH  to 
produce  methanol  (CH3OH)  and  CH3O.  Importantly,  a  recombina¬ 
tion  route  is  observed  to  form  C(2)  hydrocarbons,  which  is  a  com¬ 
mon  feature  of  the  low-temperature  oxidations.  In  the  path  of 
recombination  channel,  CH3  radicals  are  recombined  to  form 
ethane  (C2H6),  which  is  ultimately  converted  to  CO  and  CH2 
through  a  series  of  C2HX  hydrocarbons. 

CH4  reaction  pathways  for  cases  1  -3  appear  to  be  very  similar 
although  the  combustion  regime  has  been  changed  when  T\n  varies 
from  700  K  to  2700  K.  The  reason  behind  is  that  the  WSR  working 
temperature  TWsr  changes  very  little,  i.e.,  from  3090  I<  to  3338  K,  as 
Tin  is  increased.  Therefore,  according  to  ref.  18,  the  chemical  path¬ 
ways  for  the  three  cases  all  follow  the  high  temperature  reaction 
pattern.  In  contrast,  for  cases  5-6,  the  combustion  regime  changes 
from  quasi-MILD  to  MILD  when  T[n  varies  from  700  K  to  2700  K. 
Obviously,  the  CH4  reaction  pathways  differ  greatly.  This  should 
result  from  the  fact  that  Twsr  increases  significantly  from  1312  K  to 
2867  I<  as  T[n  increases  from  700  I<  to  2700  K.  Consequently,  the  low 
temperature  chemical  pathway  pattern  for  case  4  is  changed  to  the 
high  temperature  one  for  case  6  [18  .  It  is  also  interesting  to  note 
that  the  CH4  oxidation  pathways  are  similar  for  cases  1-3  and  f 
because  TWsr  ~  3000  K  for  all  these  cases.  From  the  above  dis¬ 
cussion,  although  the  investigated  cases  are  under  different  com¬ 
bustion  regimes,  the  reaction  pathways  are  found  to  mainly  depend 
on  Twsr-  In  other  words,  the  reacting  temperature  has  a  signifi¬ 
cantly  stronger  influence  on  combustion  pathways  than  does  the 
combustion  regime. 

7.  Conclusions 

The  present  study  has  characterized  and  also  classified  the  WSR 
combustion  of  G-L4/O2/N2  based  on  the  kinetic  calculations  over 
wide  ranges  of  the  inlet  mixture  temperature  (Tin),  nitrogen  con¬ 
centration  for  dilution  (XjJ,  )  and  equivalence  ratio  (<p).  The  S-curve 
is  utilized  to  identify  the  extinction  and  self-ignition  temperatures 
(Tex  and  Ts j)  of  the  CH4/02/N2  mixture,  from  which  the  mapping  of 


the  various  combustion  regimes  is  obtained.  On  the  basis  of  the 
results  presented  in  Sections  5  and  6,  several  conclusions  are  made 
below. 

1.  The  Tex  and  TS[  can  be  changed  by  varying  either  the  dilution 
degree  or  the  residence  time  in  the  reactor  or  by  both.  As  the 
dilution  is  increased  and/or  the  residence  time  is  decreased, 
both  Tex  and  TS[  increase,  but  their  difference  narrows  rapidly. 

2.  A  three-dimensional  Tin—X^  —<p  plot  is  produced.  From  it,  the 
combustion  regime  at  any  given  operative  condition  can  be 
easily  identified,  i.e.,  the  traditional  combustion  (TC),  high 
temperature  combustion  (HTC),  and  flameless  combustion 
(FLC),  where  FLC  may  be  further  divided  into  the  three  sub¬ 
zones:  MILD,  MILD-like  and  quasi-MILD. 

3.  The  classification  of  the  combustion  regime  appears  to  be 
generic  and  depend  little  on  the  combustor  configuration 
because  similar  observations  can  be  made  from  the  present  WSR 
combustion  and  diffusion  flames  under  different  configurations 
[12,13,15]. 

4.  Both  the  overall  mechanism  and  elementary  chemical  pathways 
of  CH4  oxidation  depend  mainly  on  TWsr  and  cp}  and  weakly  on 
the  combustion  regime.  This  conclusion  derives  from  analyzing 
the  fuel  conversion  to  the  composition  of  the  WSR  products 
under  various  initial  conditions. 
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